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Abstract. Pulsars are arguably the only astrophysical sources whose emission spans the en- 
tire electromagnetic spectrum, from decameter radio wavelengths to TeV energies. The LOw 
Frequency ARray (LOFAR) offers the unique possibility to study pulsars over a huge frac- 
tional bandwidth in the bottom 4 octaves of the radio window, from 15-240 MHz. Here we 
present a LOFAR study of pulsar single pulses, focussing specifically on the bright nearby pul- 
sar B0809+74. We show that the spectral width of bright low-frequency pulses can be as narrow 
as 1 MHz and scales with increasing frequency as A///c 0.15, at least in the case of the 
PSR B0809+74. This appears to be intrinsic to the pulsar, as opposed to being due to propaga- 
tion effects. If so, this behavior is consistent with predictions by the strong plasma turbulence 
model of pulsar radio emission. We also present other observed properties of the single pulses 
and discuss their relation to other single-pulse phenomena like giant pulses. 
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1. Introduction 

In recent years a number of studies of strong individual pulses from bright nearby 
pulsars were performed, including some at low radio frequencies (see, e.g., Kuzmin 2006 
and references therein; Weltevrede et al. 2006; Karuppusamy et al. 2011). Higher radio 
frequencies are more favorable for pulsar studies because of several factors, such as intra- 
channel dispersion, scattering, the Galactic synchrotron background, and the ionosphere. 
Nonetheless, the very low- frequency range is interesting for pulsar studies because of 
pulse-profile evolution, which can become quite dramatic towards the lowest frequencies 
(Hassall et al. 2012). Next, some phenomena are potentially not seen at higher frequencies 
(or become less pronounced there), in particular the anomalously bright individual pulses 
from five nearby pulsars including PSR B0809+74, reported by Ulyanov et al. (2006) 
based on observations with the UTR-2 radio telescope. It is not clear whether these 
pulses are attributed to giant pulses or other single-pulse phenomena observed from 
other pulsars at higher frequencies. In these proceedings we report the first broadband 
low-frequency single-pulse study of the PSR B0809+74 using the LOFAR telescope. 

2. Observations and Data processing 

Observations were carried out using the Low Band Antennas (LB As) in the LBA_OUTER 
array mode, where the 48 outermost dipoles in each LBA field are used. The instanta- 
neous bandwidth of 46.875 MHz centered at a frequency of about 38 MHz was split into 
240 subbands of 32 channels each and originally sampled every 491.52 /is. Each obser- 
vation was the coherent sum of the six "Superterp" stations, CS002-CS007. A detailed 
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Figure 1. Examples of characteristic spectra and profiles of individual pulses of the pulsars 
B0809+74. The number in the top-right of each panel shows the pulse number from the beginning 
of the observation. There are 438 bins in the period, and each spectrum is comprised of 480 
channels of 97.656 kHz each. 



description of LOFAR's pulsar observing modes and the online data reduction pipeline 
is given in Stappers et al. (2011). 

The data were converted from the LOFAR beam-formed format to the PSRFITS data 
format, corrected for the varying gain over the LB A band, and processed with DSPSfQ 
to form both an integrated profile as well as single-pulse integrations. We used pdmp 
from PSRCHIVeQ on a sample of strong pulses to determine the best DM for the time of 
our observations: 5.752 pc cm~^, with the DM jitter between the pulses being on the 
order of 0.002 pc cm~^. Many of the pulses showed very narrow-band structure in their 
spectra (see Figure [TJ and details below). Using a 2D time- frequency search technique 
(Kondratiev et al., in prep) we next analysed the spectra of the complete sample of single 
pulses. 



3. Results 

The complete results of our study will be presented in Kondratiev et al. (in prep); here 
we briefly present a few highlights. 

Figure [T] shows some typical pulses from the PSR B0809+74. In a given pulse, emis- 
sion can occur in either the leading or trailing component; in both components at the 
same time; or sometimes even in three components corresponding to three subpulse-drift 
bands (see, e.g., pulse #941 on Figured]). The frequency structure of the pulse spec- 
tra can be also quite different. For instance, pulse #941 shows a very narrow emission 
patch at 19 MHz in the second component, but the strongest emission from other two 
components occurs at the higher frequencies. Other examples (e.g. pulse #2003) show 
very broadband structure comprised of a number of smaller, individual patches with the 
maximum tending to be at lower frequencies (< 30 MHz). On average the pulse spectra 
of the PSR B0809+74 are quite narrow and constitute only half of the band. 

Figure [2] presents the pulse energy distribution for the pulsar B0809+74. The energy 
is calculated for the emission patches of every pulse component in the profile rather than 
in the entire band due to the presence of narrow-band pulses. Distributions are skewed 
towards higher energies due to the search technique as we searched for highest spectral 
fiux density. It is clear that distributions are very similar at low energies with slightly 
less number of low-energy patches in the ON-pulse window. A few positive outliers for 
the OFF-pulse histogram indicates the presence of RFI but their number is insignificant. 
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Figure 2. Pulse energy distribution for B0809+74. Energies are for the emission patches over 
patch's frequency width normahzed by the energy of the average profile over the entire band. 
Distributions in green and yellow are for ON-pulse and OFF-pulse phase windows, respectively, 
using the same search technique with the threshold of Icr. The distributions in light blue are for 
the strong patches with the spectral flux density over 5cr. 



We also show the ON-pulse histogram for the strong patches with the spectral peak flux 
density > 5a. The roll-off of the histogram at low energies represents the selection bias 
for patches with peak spectral fluxes close to the detection threshold. The apparent larger 
fraction of stronger patches than la-patches for the same energy, is caused by the fact 
that overall number of found patches is much larger for the la than the 5<j threshold. 

We fit a normal, log-normal and power-law distribution with low-energy cutoff. For our 
fits we excluded the very high energy tail of the distributions that have small statistics. 
The lognormal distribution provides the best goodness of fit. However, fitting the power- 
law distribution to the high-energy tail for energies > 2A{E) provides better results. In 
general, the significance of the fit is not high enough for both pulsars to reject neither 
lognormal nor power-law distributions, and the better statistics is needed. 

The narrow-band emission of the PSR B0809+74 

The narrow-band spectra of some of the pulses from the PSR B0809+74 are quite unique 
and have never been observed at higher frequencies. We found it very unlikely that 
this narrow-band frequency pattern is caused by either ISM scintillations or ionosphere. 
For the former estimates of the decorrelation bandwidth give values of < 2 kHz, so all 
scintillations should be completely averaged out in our data. For the ionosphere the rate 
of change in its state is long (mins-hours) compared to time scales of about a second for 
the use case when narrow-band pulses follow each other. 

Figure [3] shows the dependence of the patch's frequency width on their central fre- 
quency for the pulsar B0809+74 for the strong patches with the energy E > 2{E). Though 
data points are somewhat scattered, they qualitatively lie on the line, l^f j fc ^ 0.15. This 
agrees very well with the prediction from the strong plasma turbulence (SPT) model 
(Weatherall 1998). The SPT model of pulsar radio emission predicts narrow-band radia- 
tion with A/// ~ 0.1 — 0.2. If indeed true, this would also provide a direct link between 
giant pulses and anomalously intensive pulses observed at low frequencies. 



4 



V. Kondratiev et al. 



Figure 3. Dependence of the 
patch spectral width on their 
emission frequency (central fre- 
quency of the patch) for the pul- 
sar B0809+74 for the ah strong 
emission patches with the en- 
ergy E > 2{E) (green circles). 
The orange line, A///c - 0.15, 
shows the best least-square fit. 
The dotted lines show the same 
dependency but for the slope of 
0.1 and 0.2. 
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4. Summary 

We observe occasional bright narrow-band pulses from the PSR B0809+74 at frequen- 
cies 15-62 MHz. Their spectra can be as narrow as 1 MHz and tend to have a width of 
3 MHz. 

We identified pulse sequences from the pulsar B0809+74 where narrow patch of emis- 
sion is drifting up in frequency from pulse to pulse. We see evidence for similar frequency 
drift for at least one other pulsar. 

The origin of these narrow-band pulses is likely to be pulsar-intrinsic rather than due 
to propagation effects in ISM or ionosphere. 

At the moment, the observed pulse properties of the PSR B0809+74 do not allow us 
to relate low- frequency bright narrow-band pulses to "spiky emission" , or to either giant 
pulses or regular emission. 

For PSR B0809+74, the spectral width of the strong emission patches scales with in- 
creasing frequency as A///c ~ 0.15, qualitatively agreeing with the prediction of the SPT 
model. This supports the relation between bright narrow-band pulses at low frequencies 
and giant pulses. 
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